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Abstract. We demonstrate that the question whether or not a 
given topological ramified covering map of the 2-sphere is Thurston 
equivalent to a rational map is algorithmically decidable. 

1. Introduction 

This paper solves a long-standing open problem in one-dimensional 
Complex Dynamics: we show that Thurston's equivalence to a post- 
critically finite rational map is algorithmically decidable. Thurston's 
theorem [3] is central to the subject. In the case when a rational map- 
ping exists, it is essentially unique, and the proof of the theorem [3] 
supplies an iterative algorithm for computing its coefficients. When, 
for instance, the rational mapping is a quadratic polynomial, the Spi- 
der algorithm of Hubbard and Schleicher [5J, computes the coefficients 
starting from a convenient combinatorial description of the branched 
covering. However, the algorithm will go astray if the branched cov- 
ering data cannot be realized by a polynomial. Thus the question we 
answer in this work is both natural and important. 

Here is the Equivalence problem we consider in this note: 

Problem: Equivalence to a Rational Map. Given a piecewise 
linear post-critically finite ramified covering map f from C to itself, 
determine whether or not it is Thurston equivalent to a rational map. 
If it is equivalent, give an algorithm to compute the coefficients of the 
corresponding rational map (defined up to a conjugacy with a Mobius 
map). 

Our main result is: 

Theorem 1.1. The problem of equivalence to a rational map is algo- 
rithmically solvable. 



Date: May 9, 2010. 



1 



2 sylvain bonnot, mark braverman, michael yampolsky 

2. Thurston mappings 

In this section we recall the basic setting of Thurston's characteriza- 
tion of rational functions. 

2.1. Ramified covering maps. Let / : S* 2 — > S 2 be an orientation- 
preserving branched covering map of the two-sphere. We define the 
postcritical set Pf by 

Pf ■= U f on ^ 

n>0 

where Qf is the set of critical points of /. When the postcritical set Pf 
is finite we say that / is a Thurston mapping. 

Thurston equivalence. Two Thurston maps / and g are Thurston 
equivalent if there are homeomorphisms <po, <f)\ : S 2 — > S 2 such that 

(1) the maps <f> , <pi coincide on Pf, send Pf to P g and are isotopic 
rel Pf- 



(2) the diagram 



£2 S 2 
f 



commutes. 

Orbifold of a Thurston map. Given a Thurston map / : S 2 — >■ S 2 , 

we define a function Nf : S 2 — > N U oo as follows: 

Nf(x) = I 00 if x is in a cycle containing a critical point, 

lcm deg y (/ ofc ) otherwise. 

The pair (S 2 ,Nf) is called the orbifold of f. The signature of the 
orbifold (S 2 ,N f ) is the set {iV/(x) for x such that 1 < N f (x) < oo}. 
The Euler characteristic of the orbifold is given by 



X6P / 

One can prove that xl^j ^/) < 0- I n the case where x(S 2 , Nf) < 0, 
we say that the orbifold is hyperbolic. Observe that most orbifolds are 
hyperbolic: indeed, as soon as the cardinality \Pf\ > 4, the orbifold is 
hyperbolic. 
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Thurston linear transformation. We recall that a simple closed 
curve 7 C S 2 — Pf is essential if it does not bound a disk, is non- 
peripheral if it does not bound a punctured disk. 

Definition 2.1. A multicurve F on (S 2 ,Pf) is a set of disjoint, non- 
homotopic, essential, nonperipheral simple closed curves on S 2 — Pf. 
A multicurve T is f-stable if for every curve 7 G T, each component a 
°f f l {l) i> s either trivial (meaning inessential or peripheral) or homo- 
topic rel Pf to an element ofT. 

To any /-stable multicurve is associated its Thurston linear trans- 
formation /r : M r — > M r , best described by the following transition 
matrix 



where the sum is taken over all the components a of / _1 (5) which 
are isotopic rel Pf to 7. Since this matrix has nonnegative entries, 
it has a leading eigenvalue A(r) that is real and nonnegative (by the 
Perron- Frobenius theorem) . 

We can now state Thurston's theorem: 

Thurston Theorem. Let f : S 2 — >• S 2 be a Thurston map with hy- 
perbolic orbifold. Then f is Thurston equivalent to a rational function 
g if and only if X(T) < 1 for every f-stable multicurve T. The rational 
function g is unique up to conjugation with an automorphism 0/P 1 . 

When a stable multicurve T has a leading eigenvalue A(r) > 1, we 
call it a Thurston obstruction. 

Several examples. Let us first give an example of a quadratic rational 
map / with an /-stable multicurve T. The map is given by the formula: 



The picture of its Julia set is seen in Figure [fl it is popularized as the 
cover art of the Stony Brook preprint series. The map / is known as the 
mating of two quadratic Julia sets: Douady's rabbit and the basilica 



The two critical points of / are ai = and b\ = 00. Both of them 
are periodic: 




/(*) 



z 2 + c 



with c 



1 + iVS 
2 



(see e.g. 



ai = t— > a 2 1— > a 3 1— > a\ 



b\ = 00 1— > b 2 = 1 1— > b\. 
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Figure 1. A stable multicurve of a quadratic rational map. 



Our stable multicurve T consists of a single simple closed curve 7 which 
separates a^'s from b^s. It is easy to see that the corresponding tran- 
sition matrix consists of a single entry 1/2. Thus, A(r) = 1/2. 

To give an example of a Thurston obstruction, we will need to work 
a little harder. We again use the procedure known as mating. Let us 
again start with Douady's rabbit polynomial, f c {z) = z 2 + c which is 
the unique quadratic polynomial with Imc > such that the critical 
point is peridic with period 3. Thus, the postcritical set of f c is {wo = 
0, W\ = c, W2 = c 2 + c, W3 = 00}. Consider also the complex conjugate, 
the polynomial f 5 whose postcritical set we denote {w' = 0,w[ = 
c,w' 2 = c 2 + c,w' 3 = 00}. The formal mating of these two polynomials 
is the branched covering mapping of S 2 which is obtained as follows. 
We first compactify the complex plane by adjoining a circle of directions 
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Figure 2. An example of Thurston obstruction. 

at infinity, {oo-e 27Tie \ 9 e K/Z}. We denote such compactification with 
the natural topology by ©. Let us now glue two copies © l5 © 2 along 
the circles at infinity using the equivalence relation ~oo given by 



oo ■ e 2nie e ©x ~oo oo ■ e- 27Ti0 G © 2 . 



Evidently, 



S = ©j U © 2 / -oo- 5 



2 



The formal mating of / c and / g is the well-defined branched covering 
map F of the 2-sphere S which is given by f c on © 1 and f 5 on © 2 . By 
construction, this map has an invariant equator (the two circles at oo 
glued together), and its postcritical set is the union 

{wo,w u w 2 } U {w^w'^w'z}. 

An obstruction for this mapping is given by a multicurve T consisting 
of three loops ji separating Wi, w[ from the rest of the postcritical set 
(see Figure [2]). It is easy to see that Y is an F-stable multicurve, with 
the associated transition matrix 




so that A(r) = I. 
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2.2. A piecewise-linear Thurston mapping. For the purposes of 
an algorithmic analysis, we will require a finite description of a branched 
covering / : S 2 — » S 2 . 

Since we will work mainly in the piecewise linear category, it is con- 
venient to recall here some definitions. 

Simplicial complexes. Following [17] (chapter 3.2 and 3.9) we call a 
simplicial complex any locally finite collection £ of simplices satisfying 
the following two conditions: 

• a face of a simplex in E is also in S, and 

• the intersection of any two simplices in £ is either empty or a 
face of both. 

The union of all simplices in £ is the polyhedron of £ (written |£|). 
Piecewise linear maps. A map / : M — > N from a subset of an 
affine space into another affine space is piecewise linear (PL) if it is 
the restriction of a simplicial map defined on the polyhedron of some 
simplicial complex. 

We also define piecewise linear (PL) manifolds as manifolds hav- 
ing an atlas where the transition maps between overlapping charts are 
piecewise linear homeomorphisms between open subsets of M. n . It is well 
known that any piecewise linear manifold has a triangulation: there is a 
simplicial complex £ together with a homeomorphism |£| — > X which 
is assumed a PL map (see [IT], proof of theorem 3.10.2). 

One example of such a manifold is the standard piecewise linear 
(PL) 2-sphere, which is nicely described in [17] as follows: pick any 
convex 3-dimensional polyhedron K C M 3 , and consider the charts 
corresponding to all the possible orthogonal projections of the bound- 
ary (topological) sphere OK onto hyperplanes in M 3 . The manifold 
thus obtained is the standard piecewise linear 2-sphere. One can prove 
that another choice of polyhedron would lead to an isomorphic object 
(see exercise 3.9.5 in [T7]). 

It is known that in dimension three or lower, every topological mani- 
fold has a PL structure, and any two such structures are PL equivalent 
(in dimension 2, see [H], for the dimension 3 consult [1]). 
Piecewise linear branched covers. We begin by formulating the 
following proposition which describes how to lift a triangulation by a 
PL branched cover (see [3], section 6.5.4): 

Proposition 2.1 (Lifting a triangulation). Let B be a compact 
topological surface, it : X — >■ B a finite ramified cover of B. Let A be 
the set of branch points ofn, and let T be a triangulation of B such that 
A is a subset of vertices ofT (A C Kq(T) in the established notation). 
Then there exists a triangulation T' of X, unique up to a bijective 



THURSTON EQUIVALENCE IS DECIDABLE 



7 



change of indices, so that the branched covering map ir : X — > B sends 
vertices to vertices, edges to edges and faces to faces. Moreover, if 
X = B is a standard PL 2-sphere and tt is PL, then T' can be produced 
constructively given a description of T ■ 

We consider PL maps of the standard PL 2-sphere which are topo- 
logical branched coverings with a finite number of branch points. We 
call such a map a piecewise linear Thurston mapping. 

Remark 2.1 . Note that any such covering may be realized as a piecewise- 
linear branched covering map of a triangulation of C with rational ver- 
tices. An algorithmic description of a PL branched covering could thus 
either be given by the combinatorial data describing the simplicial map, 
or as a collection of affine maps of triangles in C with rational vertices. 
We will alternate between these descriptions as convenient. 

We note: 

Proposition 2.2. Every Thurston mapping f is Thurston equivalent 
to a PL Thurston mapping. 

Before proving the above Proposition, let us formulate a basic topo- 
logical fact, known as Alexander's trick: 

Alexander's trick. Two homeomorphisms of the closed n- dimensional 
ball, which are isotopic on the boundary, are isotopic. 

Proof of Proposition \2.2\ We may start with a triangulation 71 of S 2 
whose vertices include the postcritical set Pf. Refining the triangu- 
lation 7i to 72, if necessary, we isotope / to a map which leaves the 
vertices and the edges of 72 invariant. Finally, every topological map 
from a triangle to a triangle can be isotoped into a simplicial map us- 
ing Alexander's trick. We can thus further isotope our map to a PL 
Thurston mapping with triangulation 72- D 

3. Outline of the proof of Theorem 11.11 

The proof will rely on a construction of two explicit algorithms, A\ 
and A 2 , which, given a postcritically finite piecewise linear branched 
covering / : S 2 — >■ S 2 with a hyperbolic orbifold, perform the following 

tasks: 

A\ If / has a Thurston obstruction, the algorithm A\ will terminate 
and output the obstruction. It will not terminate otherwise. 

A2 If / is Thurston equivalent to a rational mapping R, then the 
algorithm A2 will terminate. It will identify the rational map- 
ping by outputting a ball in an appropriate parameter space 
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of rational maps which isolates the rational mapping R from 
postcritically finite mappings of the same degree and with the 
same size of the postcritical set. If / is not equivalent to any 
rational mapping, then the algorithm A 2 will not terminate. 
We further will use a polynomial root-finding algorithm A% which 

finds an isolated root G M. m of a system of polynomial equations 

{Pi(x) = 0}. 

A3 the input of the algorithm is: a system of polynomial equations 
{Pi(x) = 0} for x G C m (the coefficients of the polynomials 
Pi are either given through an oracle, or computed with an 
arbitrarily high precision via a given algorithm); a rational ball 
B(w,r) C C m which contains 2* and such that B(w,2r) does 
not contain any other roots; a natural number n. The output 
is d n G C m with the property — d n || < 2~ n . 

Proof of Theorem \l.l\ assuming the existence of A\ and A 2 . Given a post- 
critically finite piecewise linear map / with a hyperbolic orbifold, we 
will run the two algorithms A\ and A 2 in parallel. One and only one of 
them will terminate. If it is A%, then we conclude that / is not equiv- 
alent to any rational map. If it is A 2 then we know that a Thurston 
equivalent rational map R exists, and we are given an isolating neigh- 
borhood for it in the parameter space. The root-finding algorithm 
A 3 can then be employed to find the coefficients of R with any given 
precision. □ 

4. Some topological preliminaries 

Simple closed curves. Recall that two curves are in a minimal posi- 
tion if they realize the minimal number of intersections in their homo- 
topy classes. 

Lemma 4.1 (The Bigon Criterion). Two transverse simple closed 
curves on a surface S are in a minimal position if and only if the two 
arcs between any pair of intersection points do not bound an embedded 
disk in S. 

Let us also formulate an elementary fact: 

Lemma 4.2. Two simple closed curves on a surface S are homotopic 
if and only if they can be isotoped to boundary curves of an annulus. 

We now prove: 

Proposition 4.3. There exists an algorithm to check whether two sim- 
ple closed polygonal curves on a triangulated surface S are homotopic. 
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Proof. The algorithm works as follows: 

(I) If necessary, isotope the curves so that all the intersections are 
transverse. 

(II) While there exists a pair of intersection points which bounds a 
disk do: 

push one of the curves through the disk to remove the two 
intersection points, end do 

(III) Does there exist an intersection point? If yes, output the 
curves are not homotopic and halt. If no, proceed to step 
(IV). 

(IV) Do the two curves bound an annulus? If no, output the curves 
are not homotopic and halt. If yes, output the curves are 
homotopic and halt. 

To verify the algorithm, we note that the Bigon Criterion implies that 
step (II) can be performed until the curves are in a minimal position. 
The correctness of the algorithm now follows by Lemma 14.21 

□ 

Maps isotopic to the identity. The following theorem of Ladegail- 
lerie j8] will be useful to us in what follows: 

Theorem 4.4. Let K be a compact topological 1-complex, X an ori- 
ented compact surface with boundary, z ,ii two embeddings of K into 
the interior of X. There is an equivalence between the two following 
properties: 

(1) «o andii are isotopic by an ambient isotopy of X (fixed on dX) 

(2) i and %\ are homotopic and there is an orientation preserving 
homeomorphism h : X — >■ X such that h o i Q = i 1 . 

We formulate the following corollary: 

Proposition 4.5. There exists an algorithm A which does the follow- 
ing. Given a triangulated sphere with a finite number of punctures 
S = S 2 — Z and a triangulated homeomorphism h : S — >■ S , the algo- 
rithm identifies whether h is isotopic to the identity. 

Proof. Let x ^ Z be a vertex in the triangulation T. Consider a 
collection of closed loops ji in dT passing through the basepoint x such 
that {7j} forms a basis of iri(S) (refine the triangulation, if necessary). 
By Theorem 14.41 it is sufficient to verify that h{pfi) is homotopic to 7, 
for all i. Indeed, this is equivalent to the existence of a global isotopy 
of S which moves hfy) to 74. By the Alexander's trick, the latter 
statement means that h is isotopic to the identity. 

□ 
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Figure 3. Dehn twist of a cylinder about an equatorial 
curve 7: a vertical line I before (left) and after (right) 
the twist is shown. 

Dehn twists. Recall the definition of a Dehn twist. Let 7 be a simple 
closed curve on a surface S, and let A be a tubular neighborhood of 7. 
Choose a homeomorphism 

h : S 1 x [0, 1] A, 

which endows the annulus A with a coordinate system (6, r) where 
9 G R/mod27rZ is the angular coordinate in S 1 , and r e [0,1]. A 
Dehn twist about 7 is the homeomorphism 

/ : S 

which is identical outside A, and is given by 

/ : h(6,r) ^ h(6 + 27ir,r). 

Mapping class group. Since Thurston equivalence involves isotopies 
preserving pointwise the points of P/,we are led to consider the pure 
mapping class group PMod(S' 2 — Pf). It is the group of homeomor- 
phisms of 5* 2 — Pf fixing Pf pointwise, modulo isotopies fixing Pf 
pointwise. The mapping class group acts on the set of isotopy classes 
of simple closed curves. 
We use the following fact: 

Proposition 4.6. The group PMod(S 2 — Pf) is generated by a finite 
number of explicit Dehn twists. 

The finiteness of the number of generating twists is a classical result 
of Dehn; Lickorish [9] has made the construction explicit. See, for 
example, [5] for an exposition. 
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5. Algorithm Ai: detecting an obstruction 

Enumeration of the multicurves. We first prove the following propo- 
sition: 

Proposition 5.1. There exists an algorithm A which enumerates all 
non-peripheral multicurves on S 2 — Pf. 

For ease of reference let us state the following elementary fact: 

Proposition 5.2. Let S* denote S 2 with a finite number of punctures. 
Consider two simple closed curves 71 and 72 in S* . Assume that a 
component of S* —71 contains the same number of punctures as some 
component of S* — 72. Then there exists a self-homeomorphism of S* 
which sends 71 to 72. 

We fix a finite collection of simple closed curves c\, . . . , cm so that the 
Dehn twists 7\, . . . , Tm around those curves generate the mapping class 
group PMod(5' 2 — Pf). This construction can be performed algorith- 
mically by Proposition 14.61 We further refine our initial triangulation 
of the sphere so that these Dehn twists can be considered as piecewise 
linear maps relatively to the refined triangulation. 

For every set of j punctures with j e {2, . . . , \Pf \ — 2} we choose one 
polygonal simple closed curve which separates them from the rest of Pf. 
Denote these curves 7^ (a simple count shows that k = 1 . . . 2l p ^~ 1 — 
l^/l-l- 

Proof of Proposition I5.il To enumerate all multicurves, we proceed in- 
ductively as follows. At step 0, our collection of multicurves consists 
of all finite subsets of the set {7^}. 

At step iV G N, we generate all the possible images of the curves 7^ 
by reduced words in the Tj,Tj _1 of length less than N. Using Propo- 
sition I4.3[ we remove all duplications from this finite collection. All 
inessential or peripheral curves are likewise removed. 

We then consider all finite subsets of this collection. Using Proposi- 
tion [O] again, we remove all subsets which have been generated previ- 
ously, at steps 0, 1, . . . , N — 1. 

By Proposition 15.21 and Proposition I4.6[ every multicurve is thus 
generated. □ 

Construction of the algorithm A±. Denote A(n) the algorithm of 
Proposition l5.11 which generates the exhaustive sequence of multicurves 
r„. Set n — 1. 

(I) use Proposition 14.31 to check whether T n is invariant. If not, 
proceed to step (V). 
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(II) Compute the transition matrix M^s of the associated Thurston 
linear transformation /p n . Denote P n (X) the characteristic poly- 
nomial of M^s- 

(III) Is 1 a root of P n l If yes, go to step (VI). 

(IV) For 1 < j < n do 

— Use A3 to query whether Pj has a root A with B(X, 2~ 31 ) C 
[1 + 2~\ 00). If yes, go to step (VI). 
end do 

(V) ra-> n + 1. Return to step (I). 

(VI) Return there exists a Thurston obstruction and halt. 



6. Algorithm A 2 \ finding an equivalent rational map 

6.1. Moduli space of rational maps. Let Rat^ denote the space 
of all holomorphic maps of degree d > 2 to itself. These maps can be 
written as fractions where the polynomials p, q are relatively prime 

and d = max(deg p, deg q). It can be shown that Rat^ is a connected 
complex-analytic manifold of dimension 2d+l. Denoting Res(P, Q) the 
resultant of P and Q, one can represent Rat^ as the open set F 2d+1 /V, 
where V = {P, Q : Res(P, Q) = 0}. (See for example [16], page 169). 

Since we are interested in equivalence classes of rational maps under 
conjugation by Mobius maps, we are led to consider the moduli space 
M,d = Ratd/PSL 2 (C) (observe that it is the same as Ratd/PGL 2 (C)). 

We note the following easy fact: 

Proposition 6.1. The moduli space Aid has the structure of a complex 
orbifold of dimension 2d — 2. 

Proof. The stabilizer of a rational map f ^ Id under the action of 
PSL 2 (C) is the subgroup S(f) consisting of Mobius maps which com- 
mute with /. There exists n 6 N such that the set P n consisting of 
periodic points of / with periods less than n has at least three points. 
Since every M e <S(f) must permute the points in P n , the stabilizer 
S(f) is necessarily finite. □ 

As an example of a rational map with a non-trivial stabilizer, con- 
sider f(z) = z d for d > 2, in which case, 

S(f) = {z^ Xz\ where X^ 1 = 1}. 

We will now require a more computation-friendly description of Aid- 
There are several similar approaches to this in the existing literature; 
we use the work [2]. As a first step we note the following standard fact: 



THURSTON EQUIVALENCE IS DECIDABLE 



13 



Proposition 6.2. Suppose R is not conjugate to a map of the form 
z i — y z ±d . Then the union of the critical and the postcritical sets CfUPf 
contains at least three points. 

Mapping scheme. A mapping scheme of degree d is a triple (N, r, u), 
where N e N and N > 3; r is a dynamics function 

{l,...,jV}-^{l,...,JV}, 

and u is a local degree function 

{1,...,N}^N. 

For a postcritically finite branched covering map / : S 2 — > S 2 denote 
Zf is the union C/UP/ of the critical set and the postcritical set of 
/. We say that / realizes X = (N, r, u) if we can choose a bijection 
V> : {1, • • • , N} Z f such that: 
. rl>(r(x)) = /(^(ar)); 

• the local degree of / at ip(x) is equal to ui{x). 

Following [2j, a normalization of a mapping scheme X = (N,t,oj) is 
an injection 

a: {0,l,oo} ->{1,...,JV}. 
We will denote a pair (X, a) by X a , and refer to it as a marked mapping 
scheme. 

We say that a rational map / : C — > C realizes a marked map- 
ping scheme X a if the bijection ^ as defined above has the additional 
property: 

• {0, 1, oo} C Zf and ip^ii) = a(i). 

The set of all rational realizations of a specific marked mapping scheme 
X* = X a will be written as Rat x (X*). 

Viewing the points in {1, . . . , N} as vertices of a weighted directed 
graph with arrows connecting x with t(x) having weights uj(x), we 
define a signature Afx as follows. The signature is the set 

Af x = {N(x), xer({l,...,N})}, 

where N(x) is defined as the least common multiple, over all directed 
paths of any length joining y to x, of the product of the weights of 
edges along this path. The significance of this definition lies in the 
following: 

Proposition 6.3 Q3J). If f is a rational mapping which realizes the 
mapping scheme X, then Afx is the signature of the orbifold of f. 

We can now state: 
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Theorem 6.4 ([2]). Given a mapping scheme X whose signature is not 
(2, 2, 2, 2) and a normalized marking scheme X* = (X, r, u), there is an 
injection i : Rat x (X*) — > C N such that the image is a zero- dimensional 
affine variety V(X) determined by an ideal X = Xx* , where X is defined 
over the rationals. In particular, Rat x (X*) is finite. Furthermore, a 
basis for X can be algorithmically computed. 

Let us give an indication of how the injection i may be defined. 
Enumerate the elements of Zf as {z{\ so that ip(i) = Z{. Isolate the 
elements sent respectively to zero and the infinity as follows: 

FZ = {n\z n G Zf - {oo} and f(z n ) = 0}, 

and FP = {m\z m G Zf — {oo} and/(^ m ) = oo}. 
Then any normalized rational realization / can be written uniquely as 

(6 1) f(z) — a o + • • • + a r z r IlneFz( z ~ z n) 



b + b lZ + . . . + 1 • Z s UmeFp( Z ~ Z m) dm 

Thus, the coordinates dj, bj, and Zj, specify the normalized rational 
map as a point in C N . 

6.2. The algorithm A 2 . We start with a piecewise-linear Thurston 
map / : S 2 —> S 2 with triangulation T. Let X be a mapping scheme 
which is realized by /, and let a be a marking of X. Set X* = X a . Let 
N and Xx* be as in Theorem 16 .41 Compute the finite set g — (gi, ■ ■ ■ g r ) 
of polynomials with rational coefficients, which generates the ideal Ix* ■ 
Denote A the algorithm of Proposition 15.11 Let {Ri, . . . , R m } be the 
finite collection of rational maps which realize the marked mapping 
scheme X*. 

The algorithm A 2 works as follows: 

(I) use A to enumerate as D\ the representatives of the mapping 

class groups PMod(S' 2 — -Pr*); 
(II) for every 1 < k < m and every pair C N x N do 

— discretize the map {Dj) -1 o R k o Dj to a piecewise linear 
Thurston mapping Mj with triangulation 71. Note that 
at this stage the algorithm A 3 may need to be invoked to 
better estimate the coefficients of R^, 

— refine the triangulations T and 71 to obtain a triangulation 
% on which both / and Mjj are defined; 

— identify all triangulated orientation preserving homeomor- 
phisms hj of %; use the algorithm of Proposition 14.51 to list 
all Wk = hj k which are isotopic to the identity; 
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— perform a finite check to determine whether there exists a 
pair Wi, Wj such that 

Wi o Mfj = wj o /. 

If yes, go to (III). End do 
(III) output / is Thurston equivalent to a rational map, output 
the isolating neighborhood for Rk, and exit the algorithm. 

7. Concluding remarks 

Let us note an easy corollary of our main result. Consider the fol- 
lowing decidability problem: 

Problem (a): Given two piecewise linear Thurston mappings f and 
g with hyperbolic orbifolds and without Thurston obstructions, are f 
and g Thurston equivalent? 

Theorem 7.1. Problem (a) is algorithmically decidable. 

Proof. Denote Rf and R g the rational maps equivalent to / and g 
respectively. The existence of such maps is guaranteed by Thurston's 
theorem. They are defined up to a Mobius conjugacy. The algorithm 
works as follows: 

(I) Check if / and g have identical mapping schemes X = (N, r, u). 
If not, output the maps are not Thurston equivalent and 

exit. 

(II) Run the algorithm A 2 to find isolating 2 _ra -neighborhoods Uf 
and U g of Rf and R g respectively in the parameter space of 
normal forms (16. ip . 

(III) For each ordered triple of distinct natural numbers (a, b, c) be- 
tween 1 and N do 

— Normalize the mapping scheme X by a : (0,1, oo) i— > 
(a,b,c). Calculate a point of the parameter space W( a ,b,c) 
representing the normal form ( 16. ip of Rf corresponding to 
(X, a) with precision 2~ ( - n+2 \ 

— If wufrA G U g then output the maps are Thurston 
equivalent and exit. 

(IV) Output the maps are not Thurston equivalent and exit. 

□ 

Consider the following natural generalization of Problem (a), sug- 
gested to us by M. Lyubich: 

Problem (b): Given two piecewise linear Thurston mappings f and 
g with hyperbolic orbifolds, are f and g Thurston equivalent? 
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Algorithmic decidability of Problem (b) presents an interesting direc- 
tion of further study. 

We find that recognizability problems of combinatorial equivalence 
of piecewise-linear maps are analogous to recognizability problems of 
piecewise-linear manifolds, which are solvable in dimensions 3 or less, 
and become algorithmically intractable in dimensions greater than 4 
(and possibly 4 as well) - see, for example, the book of Weinberger 
[TB] • In conclusion, we speculate that natural notions of equivalence 
of maps in higher dimensions will lead to algorithmically unsolvable 
problems - and a new interplay between Dynamics and Computability. 

Appendix A. A root-finding algorithm. 

The existence of a root-finding algorithm A 3 is a classical result of 

H. Weyl [19J. Consider a system of analytic functions g = {ft (z), i = 

I, . . . ,n} defined in a box D C C n . Denote N(d,D) the number of 
common zeroes of g in D, and assume that there are no common zeroes 
on the boundary. Then N(d, D) can be computed using the following 
multidimensional residue formula (see [15] page 324): 

(ti 1)! f 1 n ^.a^ 

N (g, D ) = , n -x n ' / r-|w ^gjdgjAdgjAdgiAdgi. . . [j] ... d~g~Adg r , 
\£Ki) j dD \g\ . =1 

where [j] means that we omit the term dgj A dg~j. 

Observe that our space Rat x (X*) is indeed given as the zeroes of N 
polynomials in C^. 

Weyl's algorithm to locate all roots of g = in an isolating neigh- 
borhood G works as follows. Begin by setting j = covering G by a 
cubic grid of size 2~- ? = 1. In each cube Ck of the grid, use the residue 
formula to check whether there are any zeros in it. Since we cannot 
catch zeros on the boundary of a cube, perform the check for a cube of 
twice the size - it is guaranteed to catch any zeros in the closure of C^. 

Throw away all cubes without any zeros. Increment j t— > j + 1 and 
divide the remaining cubes Ck into cubes with side 2~- J . Repeat the 
process, until all zeros are identified with the desired precision. 
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